ON THE ISING MODEL WITH RANDOM BOUNDARY CONDITION 
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Abstract. The infinite-volume limit behavior of the 2d Ising model under pos- 
sibly strong random boundary conditions is studied. The model exhibits chaotic 
size-dependence at low temperatures and we prove that the '+' and '-' phases are 
the only almost sure limit Gibbs measures, assuming that the limit is taken along 
a sparse enough sequence of squares. In particular, we provide an argument to 
show that in a sufficiently large volume a typical spin configuration under a typ- 
ical boundary condition contains no interfaces. In order to exclude mixtures as 
possible limit points, a detailed multi-scale contour analysis is performed. 



1. Introduction 

A fundamental problem in equilibrium statistical mechanics is to determine the 
set of physically accessible thermodynamic states for models defined via a family 
of local interactions. Usually [15, 23] one interprets the extremal elements of the 
set of translationally invariant Gibbs measures as the pure thermodynamic phases 
of the model. In particular this means that one gathers all periodic or quasiperi- 
odic extremal Gibbs measures into symmetry-equivalent classes and identifies the 
latter with the pure phases. Examples are the ferromagnetic, the antiferromag- 
netic, crystalline or quasicrystalline phases exhibited by various models. In this 
approach one does not consider either interface states or mixtures as pure phases. 
The mixtures allow for a unique decomposition into the extremal measures and 
are traditionally interpreted in terms of a lack of the knowledge about the ther- 
modynamic state of the system. They can also be classified as less stable than 
the extremal measures [25, 36]. It is thought that interface states which are ex- 
tremal Gibbs measures are more stable than mixed states, but less so than pure 
phases. However, such an "intrinsic" characterization has not been developed. 
Note, moreover, that in disordered systems such as spin glasses, the stability of 
pure phases is a priori not clear and characterizing them remains an open ques- 
tion. 



1991 Mathematics Subject Classification. 82B20, 82B44, 60F05. 

Key words and phrases. Random boimdary conditions, metastates, contour models, multi-scale 
analysis, local-limit theorems. 

1 



2 



A. C. D. VAN ENTER, K. NETOCNY, AND H. G. SCHAAP 



An efficient strategy for models with a simple enough structure of low-tempera- 
ture phases is to associate these with suitable coherent boundary conditions. The 
latter are usually chosen as ground states of the model. As an example, the '+' and 

Ising phases can be obtained by fixing the constant '+', respectively the constant 
'-' configurations at the boundaries and by letting the volume tend to infinity 
This idea has been generalized to a wide class of models with both a finite and a 
'mildly' infinite number of ground states, and is usually referred to as the Pirogov- 
Sinai theory [4, 11, 49, 45, 46]. The main assumption is that the different ground 
states are separated by high enough energy barriers, which can be described in 
terms of domain walls, referred to as contours. A useful criterion to check this so- 
called Peierls condition is within the formalism of m-potentials due to Holzstynski 
and Slawny [27]. 

An alternative strategy is to employ a boundary condition that does not favor 
any of the phases. Examples are the free and periodic boundary conditions for the 
zero-field Ising model, or the periodic boundary conditions for the Potts model 
at the critical temperature. In all these cases, an infinite-volume Gibbs measure is 
obtained that is a homogenous mixture of all phases. 

Another scenario has been expected to occur for spin glasses. Namely, New- 
man and Stein have conjectured [37, 38, 39, 41, 42] that some spin glass models 
under symmetric boundary conditions exhibit non-convergence to a single ther- 
modynamic limit measure, a phenomenon called chaotic size dependence (see also 
[19, 34, 14]). In this case, both the set of limit points of the sequence of the finite- 
volume Gibbs measures and their empirical frequency along the sequence of in- 
creasing volumes are of interest, and the formalism of metastates has been devel- 
oped [39, 41, 40] to deal with these phenomena. These arguments have been made 
rigorous for a class of mean-field models [8, 30, 32, 7, 17, 43, 31], whereas no such 
results are available for short-range spin glasses. For some general background 
on spin glasses and disordered models we refer to [6, 20, 33, 47]. 

A natural toy-problem where the usual contour methods can be used in the 
regime of chaotic size-dependence is the zero field Ising model with the boundary 
condition sampled from a random distribution which is symmetric under the spin 
flip. In dimension 2 or more and at any subcritical temperature (including T = 0) 
the finite-volume Gibbs measures are expected to oscillate randomly between the 
'+' and the '-' phases, demonstrating the chaotic size dependence with exactly two 
limit points coinciding with the thermodynamic phases of the model [38]. In par- 
ticular, one does not expect either any interface (e.g. Dobrushin) Gibbs states or 
any non-trivial statistical mixtures to occur as the limit points. This problem was 
addressed in [16] where the conjecture was rigorously proven as the almost sure 
picture in the regime of the weak boundary coupling. In this regime, the bound- 
ary bonds are made sufficiently weaker w.r.t. the bulk bonds so that the interface 
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conjfigurations become damped exponentially with the size of the system, uni- 
formly for all boundary conditions. Hence, all translationally non-invariant Gibbs 
measures are forbidden as possible limit points and one only needs to prove that 
the mixtures do not appear with probability 1. 

In this paper we continue this study by removing the weakness assumption on 
the boundary bonds. To be specific, we consider the 2d Ising model with the ran- 
dom boundary condition sampled from the symmetric i.i.d. field {—1, 1}^^ and 
coupled to the system via the bulk coupling constant. The conjecture remains true 
in this case and the crucial novelty of our approach is a detailed multi-scale analy- 
sis of contour models in the regime where realizations of the boundary condition 
are allowed that violate the 'diluteness' (Peierls) condition, possibly making in- 
terfaces likely configurations. To be precise, these interfaces can have large Gibbs 
probabilities for certain boundary conditions, but we will show that such bound- 
ary conditions are sufficiently unlikely to occur for large volumes. An important 
side-result is the almost sure absence of interface configurations. This means that 
for a typical boundary condition, the probability of the set of configurations con- 
taining an interface tends to zero in the infinite-volume limit. Note that this ex- 
cludes interfaces in a stronger way than the familiar result about the absence of 
translationally non-invariant Gibbs measures in the 2d Ising model [22, 1, 26]. In- 
deed, the absence of fluctuating interfaces basically means that not only the expec- 
tations of local functions but also their space averages (e.g. the volume-averaged 
magnetization) have only two limit points, corresponding to the two Ising phases. 
Hence, we believe that our techniques allow for a natural generalization to any di- 
mension d >2. However, as already argued in [16], in dimensions d > A, the set 
{/i"*", is expected (and partially proven) to be the almost sure set of limit mea- 
sures, the limit being taken along the regular sequence of cubes. On the other 
hand, for d = 2,3 the same result can only be obtained if the limit is taken along 
a sparse enough sequence of cubes. In the latter case it remains an open problem 
to analyze the set of limit points along the regular sequence of cubes. Our con- 
jecture is that the almost sure set of limit points coincides then with the set of all 
translationally invariant Gibbs measures, i.e. including the mixtures. 

The structure of the paper is as follows. We will first introduce our notation in 
section 2, and describe our results in section 3. Then in sections 4 and 5 we will 
introduce a contour representation of the model and set up our cluster expansion 
formalism. In section 6 we first exclude the occurrence of interfaces. In the rest 
of the paper we develop a multiscale argument, providing a weak version of the 
local limit theorem to show that no mixed states can occur as limit points in the 
infinite-volume limit. Two general results, the first one on a variant of the cluster 
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expansion convergence criteria for polymer models and the second one on local 
limit upper bounds, are collected in two Appendices. 

2. Set-up 

We consider the two-dimensional square lattice T? and use the symbols a, 77, . . . 
for the maps 1? 1— >■ {—1, 1}. They are called spin configurations and the set of all 
spin configurations isQ = {— Furthermore, the symbol a a is used for the 
restriction of a spin configuration a G to the set A all? .\i A — {x), we write 
instead. The set of all restrictions of ^2 to the set A is ^a- 

A function f : Q Ris called local whenever there is a finite set D d 7? such 
that (Td — (j'd implies /(cr) = /(cr'). The smallest set with this property is called 
the dependence set of the function / and we use the symbol Vf for it. To every local 
function / we assign the supremum norm ||/|| = sup^g^ 1/(^)1- 
The spin configuration space comes equipped with the product topology, which 
is followed by the weak topology on the space M{Vt) of all probability measures 
on Q. The latter is introduced via the collection of seminorms 

\\^\\x = sup (1) 
11/11=1 

■DfCX 

upon all finite X <zZ^. Then, the weak topology is generated by the collection of 
open balls Bx{ij,) — {i/; ||z/ — < e}, e > 0, X finite, and a sequence e M(Q) 
weakly converges to fi if and only if — for all finite X c Z^. Under 

the weak topology, M{Q) is compact. 

We consider a collection of the Hamiltonians H'J^ : Q\ M. for all square vol- 
umes A = A(A^), = 1, 2, . . ., 

A{N) = {x e Z^; \\x\\ < N} = maxdxil, |a;2|} (2) 

and boundary conditions rj efl. The Hamiltonians are given by 

{x,y)cA A^^'^Lc 

where {x, y) stands for pairs of nearest neighboring sites, i.e. such that \\x — y\\i ■= 
ki — yi| + ^2 — I/2I = 1/ and A'^ = \ A. We consider the ferromagnetic case, (3 > 
0. Following a familiar framework, we introduce the finite-volume Gibbs measure 
Hl e M(Q) by 

fJ^W = ^ CXp[-HI{(Ja)] t{aj,c=vj,c} (4) 

and define the set Qp of (infinite-volume) Gibbs measures, Qp, as the weak closure of 
the convex hull over the set of all weak limit points of the sequences {i/^j^~^)n ^00, 
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T) eD.. A standard result reads that there exists (3c such that for any /3 > /3c the set of 
Gibbs measures Gfs = {a^'^ + (1 — a)ix~; < a < 1}. Here, the extremal measures 
/x^ are translation-invariant, they satisfy the symmetry relation / d/i^(cr) / (a) = 
/ dyLi~((T) /(— cr), and can be obtained as the weak limits limjv^oo A*A(Ar) V = 

3. Results 

We consider the limit behavior of the sequence of finite-volume Gibbs mea- 
sures (At^(jv)) ^eN under boundary conditions rj sampled from the i.i.d. symmetric 
random field 

P{rj, = 1} = P{ri, = = ^ (5) 

Our first result concerns the almost sure structure of the set of all limit points of 
the sequence of the finite-volume Gibbs measures, the limit being taken along a 
sparse enough sequence of squares. 

Theorem 3.1. For arbitrary uj > there is a (3i ~ such that for any (3 > (3i the set 
of all weak limit points of any sequence (/iA(fcjv))Af=i,2,.../ kN > N"^^"^, is /i"}, P-a.s. 

Remark 3.2. The above theorem does not exclude other measures as the almost sure limit 
points, provided that other (non-sparse) sequences of squares are taken instead. Actually, 
our conjecture is that, for (3 large enough, the set of all weak limit points o/(/XA(iv))iv=i,2,... 
coincides P-a.s. with Qp. On the other hand, in dimension 3, it is rather expected to 
coincide with the set of all translation-invariant Gibbs measures, and, in any dimension 
higher than 3, with the set 

Remark 3.3. A modification of the Hamiltonian (3) is obtained by re-scaling the boundary 
coupling by a factor A to get 

In this case, the claim of Theorem 3.1 for the sequence of the finite-volume Gibbs measures 

/^a''(^) = i exp[-i7^'''(aA)] ^{a^c=VAc} (7) 

was proven in [16] under the condition that \X\ is small enough (= the boundary coupling 
is sufficiently weak w.r.t. the bulk one). It was also shown that is the almost 

sure set of limit points of the sequence (/x^^j^^ atgN/ provided that the space dimension is 
at least 4. 



To reveal the nature of all possible limit points that can appear along the se- 
quence of squares A(iV), N = 1,2,..., we study the empirical frequency for the 
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jfinite-volume Gibbs states from the sequence (/x^(^))iveN to occur in a jfixed set of 
measures. More precisely, for any set B c M(Q), boundary condition 77 e and 
N — 1,2, . . ., we define 

1 ^ 

= nJ:m,,^b, (8) 

k=l 

The next theorem shows the null-recurrent character of all measures different from 
both /x"*" and We use the notation B and B^ for the weak closure and the weak 
interior of B, respectively. 

Theorem 3.4. There is /?2 such that for any f3 > (^2 any set B c MiVt), one has 

'0 ifiJ,+,i2-^B 
lim g^''' = <j i iffi^ e and fi^ ^ B (9) 



1 ifii+,ii-eB^ 



with P-probability 1. 

Both theorems follow in a straightforward way from the following key estimate 
that will be proven in the sequel of the paper. 

Proposition 3.5. Given a > 0, there is a (5q — /?o(«) such that for any (5 > (5q, e > Q 
and X c finite, 

^^m-^ p{{\\^x\^^ - ^x^x a - ^x-\\x) > e} < 00 (lO) 

Remark 3.6. The proposition claims that, for a typical 77 e the finite-volume Gibbs 
measures are expected to be near the extremal Gibbs measures /j,^. The above probability 
upper-bound of the form 0[N^2+°^ will be proven by means of a variant of the local limit 
theorem for the sum of weakly dependent random variables. Although we conjecture the 
correct asymptotics to be of order N~2, the proof of any lower bound goes beyond the 
presented technique. This is why the detailed structure of the almost sure set of the limit 
Gibbs measures is not available, except for the limits taken along sparse enough sequences 
of squares. 

Proof of Theorem 3.1. Given uj > 0, we choose an a < a;/(2(2 + u)) and define 
(3^{uj) = l3o{a). Let f3 > /3i(cj) and > iV2+^. 

First let fi ^ {/x"*", There exists a weakly open set B c M{^1) such that /x e 5 
and /i" ^ B. Choosing a finite set X (zl? and e > such that B\{\x^) r\B = %, 
Proposition 3.5 gives the bound 

= 0(A;(iV)-i+«) = 0(Ar-i+"(2+-)-f ) 
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Since P{ii\^^^^-^ e B} < oo, the set B contains P-a.s. no limit points of the 
sequence ^J^^^k^) due to the Borel-Cantelli argument. Hence, with P-probability 1, 
/X is not a limit point. 

To prove that both /i+ and /i are P-a.s. limit points, take any finite set of sites 
X and e > such that B%^{^^) fi fi^c (/" ) = 0- By the symmetry of the distri- 
bution, P{/iA(fcjv) e = -P{/^A(feiv) e B%^{ii-)] and, employing Proposi- 
tion 3.5 again, limjv P{iJ.A(k]y) e ^x(a*'^)} = |- By the Borel-Cantelli and the com- 
pactness arguments, the weak closure Bx{ij,^) contains a limit point, P-a.s. As 
/j,^ — r\x,eBx{lJ'^), the statement is proven. □ 

Proof ofTheorem 3.4. Choose ^2 — Po{oi) for an arbitrary a e (0, |) and assume 
P > p2, B e M{n). Using the notation g^'" = P{a*I(jv) ^ B} and repeating the 
reasoning in the proof of Theorem 3.1, one gets 



1 _ o{N--2+'') ^1 if e 5° and ^ B 



(12) 



and 



Hence, Ejv ]^ Var l^^^^^^g^} < oo and since the functions = 1, 2, . . . 

are independent, the result immediately follows from the strong law of large num- 
bers [13]. □ 

4. Geometrical representation of the model 

We define the dual lattice {I?)* = 1? + (1/2, 1/2). The (unordered) pairs of near- 
est neighboring sites {x, y) c 1? are called bonds and to every bond we assign a 
unique dual bond {x*,y*) = {x, y)* C (Z^)*. Given a set of dual bonds A*, we use 
the symbol \A*\ to denote the number of all dual bonds in A*. Further, with a 
slight abuse of notation, we also write x* e A* whenever there exists a dual bond 
{x*, y*) e A*, i.e. A* also stands for the corresponding set of dual sites. 

Any set A* of dual bonds is called connected whenever for any dual sites x*,y* e 
A* there exists a sequence of dual bonds {x*, x^), {xl, x^),. . . , {xl_^, y*) e A*. The 
distance d[A*, B*\ of the sets of dual bonds A*, B* is defined as the smallest in- 
teger k such that there exist x* E A*, y* G B*, and a sequence of dual bonds 
{x*,xl), {xl, X2),. . . , {xl_i, y*) C {I?)*. Similarly, a set of sites A dl? is called con- 
necteti whenever for all a;, 1/ e /I there exists a sequence of bonds {x,xi), {xi,X2), . . . , 
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Figure 1 . Pre-contours constructed via the rounding corner procedure. 

{xk-i, y) C A. Correspondingly, the distance d[A, B] of the sets A,B a I? is un- 
derstood in the sense of the || . || i-norm. 

In the sequel we assume that a volume A = A(A^) is fixed and we define the 
boundary dA as the set of all dual bonds {x, y)* such that x e A and y e A''. In 
general, dA, ^4 c A is the set of all dual bonds {x,y)*, x e A, y e A'^. For any 
subset P c dA we use the symbol P to denote the set of all sites y e A'^ such that 
there is a (unique) bond (,x, y)* E P, x e A. If P is a connected set of sites, then 
P is called a boundary interval. Obviously, any boundary interval is a connected 
set of dual bonds, however, the opposite is not true. However, any set P a dA 
has a unique decomposition into a family of (maximal) boundary intervals. Fur- 
thermore, consider all connected sets Pj of dual bonds satisfying P C Pj C dA 
which are minimal in the sense of inclusion. The smallest of these sets is called 
Con(P) (in the case of an equal size take the first one in the lexicographic order) 
and we use the shorthand |P|con = I Con(P)|. Finally, we define the comers of 
A{N) as the dual sites ^ = (-A^ - 1/2, -A^ - 1/2), ,x^2 = (N + 1/2, -A^ - 1/2), 
a;^^3 = (A^ + 1/2, A^ + 1/2), and a;^_4 = (-A^ - 1/2, A^ + 1/2). 

Pre-contours. Given a configuration a e Qa — {—1, +1}^, the dual bond {x,y)* 
to a bond {x, y) c A is called broken whenever ax ^ o-y, and the set of all the 
broken dual bonds is denoted by A\{a). In order to define a suitable decomposi- 
tion of the set Aa (cr) into components, we take advantage of a certain freedom in 
such a construction to obtain the components with suitable geometrical proper- 
ties. In this first step, we define the pre-contours as follows. Consider all maximal 
connected components of the set of dual bonds AA(cr). By the standard 'rounding- 
comer' procedure, see Figure 1, we further split them into connected (not neces- 
sarily disjoint) subsets, 7, which can be identified with (open or closed) simple 
curves. Namely, 

7 = {{4, xl), {xl,x;),..., {xl_^, xl)} keN 
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A 

Figure 2. Small boundary pre-contour. 



such that if x* = x*,i j, then {i, j} = {0, k} and 7 is closed. Otherwise, x* ^ x* 
for all i ^ j and 7 is open with Xq, xl E dA. 

These 7 are called pre-contours and we use the symbol ^'a(o-) for the set of all 
pre-contours corresponding to a; write also = {'^a(o'), cr e Ha} and use 
the symbol for the set of all pre-contours in A. Any pair of pre-contours 
7i,72 e ^A is called compatible whenever there is a configuration a e Qa such 
that 71, 72 e ^a((7). a set of pairwise compatible pre-contours is called a compati- 
ble set. Obviously, ^a is simply the collection of all compatible sets of pre-contours 
from )Ca- Intuitively, the pre-contours that are closed curves coincide with the fa- 
miliar Ising contours, whereas the pre-contours touching the boundary become 
open curves. 

Obviously, Qa i-^ ^a is a two-to-one map with the images of the configurations 
(7 and —a being identical. In order to further analyze this map, we introduce the 
concept of interior and exterior of the pre-contours briefly as follows (the details 
can be found in [5, 16]). If a G CIa is a configuration such that VAicr) = {7}, 
then there is a unique decomposition of the set A into a pair of disjoint connected 
subsets, A = Ai U A2, such that for any bond {x,y), x e Ai, y e A2, one has 
{x, y)* e 7. These are called the exterior, Ext(7), and the interior, Int(7), where the 
assignment is given by the following procedure. We distinguish three mutually 
exclusive classes of pre-contours: 

i) Bulk pre-contours. 

dA = dAi. Then, Ext (7) := Ai and Int(7) := A2, see Figure 1. 
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ExtO) 






Int(t) 




Figure 3. 


Interface. 



A 



ii) Small boundary pre-contours. 

Ai contains at least three corners of A and 9A2 7^ 0. Then, Ext(7) := Ai and 
Int(7) := A2, see Figure 2. 

iii) Interfaces. 

Both Ai and A2 contain exactly two comers of A and a) |Ai| > IA2I, or b) 
l-'^il = 1-^2 1 and xc,i e Ai. Then, Ext(7) := Ai and Int(7) := A2, see Figure 3. 

The set :— 9Int(7) is called the boundary of the pre-contour 7. 

Contours. Next, we define contours by gluing some boundary pre-contours to- 
gether via the following procedure. Any compatible pair of pre-contours 71, 72 G 

iC\is calledboundary-matching iff d'yirid'y2 ^ 0- Any compatible set of pre-contours 
such that the graph on this set obtained by connecting the pairs of boundary- 
matching pre-contours becomes connected is called a contour. In particular, ev- 
ery bulk pre-contour is boundary-matching with no other compatible pre-contour. 
Therefore, every bulk pre-contour is trivially a contour. We use the symbol VA{cr) 
for the set of all contours corresponding to a e Qa and /Ca for the set of all con- 
tours in A. Any pair of contours Fi, r2 is compatible, Fi ~ F2, whenever all pairs 
of pre-contours 71 G Fi, 72 G F2 are compatible, and we write for the set of all 
families of pairwise compatible contours in A. All the above geometrical notions 
naturally carry over to contours and we define the exterior, Ext(F) :— ri-ygr Ext(7), 
the interior, Int(F) A \ Ext(F) (in general, not a connected set anymore), the 
boundary dV := U-^er<97, and the length |F| := J^-yev bl- Similarly, if (9 G is 
a configuration of contours, let Ext(5) := nreaExt(F), lnt{d) := A \ Ext(5), and 
1^1 := Ereon 
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Eventually we arrive at the following picture. The set /Ca of contours is a union 
of three disjoint sets of contours, namely of the sets of all 

i) bulk (pre-) contours. 

ii) small boundary contours V defined by 1) 9r 0, and 2) no pre-contour 7 e F 
is an interface. 

a) simple small boundary contours: the boundary dV contains no corner, 
i.e. dV is a boundary interval. 

b) corner small boundary contours: there is exactly one comer ^ e dV. 

iii) large boundary contours V, i.e. containing at least one interface 7 e F. 

Examples of the bulk, small boundary, and large boundary contours are given in 
Figures 4-5. 

Furthermore, I'a(o') is a two- to-one map 1-^ ^a satisfying the spin-flip sym- 
metry r'A(c) = T^k{—o'). Since a takes a unique spin value in the set Ext(r'A(c"))/ 
there is a natural decomposition f^A = ^^a ^ according to this value, i.e. 



^ f|Ext(X>A((T)) = ±1} 



(14) 



As a consequence, Pa splits into a conjugated (by spin-flip symmetry) pair of 
one-to-one maps Vt^ ^ ^a- This enables us to represent the finite-volume Gibbs 
measure (4) in the form of a convex combination of two conjugated constrained 
Gibbs measures as follows: 



1 + 



^A 



(^) + 



1+ 



^A 
^A 



'''(a) 



(15) 
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Figure 5. Large boundary contour. 

where we have introduced the Gibbs measure constrained to Q J by 

= i ^M-Hl{a)] (16) 

Moreover, for any a e Q^, the Hamiltonian can be written as 

Hl{a) = Ei^\d) + 2pY.\^\ (17) 

Ted 

with d — Vj^{a), and we have introduced 

E^^\d) = -(3 Yl <^^Vy (18) 

Finally, Z^'^' is essentially the partition function of a polymer model [29], see also 
Appendix A, 

deS>A Ted 

where the polymers coincide with the contours and the polymer weights are de- 
fined by 

p±'''(r) = exp (-2/3|r|) exp {-E^^^iT) + ^^'"(0)) (20) 

By the spin-flip symmetry, we can confine ourselves to the case and use the 
shorthand notations p''(r) := ^+'"(1) = p-'-^lF), Zl := = Z^'"^, El := 

= E]^'-", and u'lia) := ly^'^^ia) = Moreover, the boundary dT of a 
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contour F has a natural decomposition into components as follows. Let cr e ilj be 
such that Vj>^{a) = {T}. Then the '±' boundary component dT^ is defined as the 
set of all dual bonds {x, y)* such that x e A, y e A^, cr^ = ±1. With this definition, 
the contour weight (20) is 

p''(r)=exp[-2/?(|r|+ ^ r^,)] (21) 

Using the representation (15) of the finite-volume Gibbs measure iJi\, the strat- 
egy of our proof consists of two main parts: 

(1) To prove that the constrained (random) Gibbs measure i>\ asymptotically 
coincides with the Ising phase, for almost all -q. 

(2) To show that a sufficiently sparse subsequence of the sequence of random 
free energy differences \ogZ\ — logZ]^^ has +00 and —00 as the only limit 
points, for almost all 77. 

Then, Proposition 3.5 follows almost immediately. Moreover, we will show that 
for a P-typical boundary condition 77 and a //^-typical configuration a e the 
corresponding set of pre-contours t>j^{a) contains no interfaces. 

Theorem 4.1. There is /S^ such that for any (5> (5^ one has 

lim AiL^^^{^A(A^)(c7) contains an interface} — (22) 

for P-a.e. rj e ^l. 

Remark 4.2. Note that the low-temperature result by Gallavotti [22], extended to all sub- 
critical temperatures in [1, 26], about the absence of translationally non-invariant Gibbs 
measures in the Id Ising model does not exclude fluctuating interfaces under a suitably 
arranged CDobrushin-like') boundary condition. On the other hand, the above theorem 
claims that a typical boundary condition gives rise to a Gibbs measure in which inter- 
faces anywhere are suppressed. We mention this side-result to demonstrate the robust- 
ness of the presented multi-scale approach and to argue that it is essentially dimension- 
independent, the d — 2 case being chosen only for simplicity. 

It is easy to realize that, for a t3^ical 77, the polymer model (19) fails the 'di- 
luteness' condition on the sufficient exponential decay of the polymer weights, 
which means one cannot directly apply the familiar formalism of cluster expan- 
sions. These violations of the diluteness condition occur locally along the bound- 
ary with low probability, and hence have typically low densities. Nevertheless, 
their presence on all scales forces a sequential, multi-scale, treatment. Multi-scale 
methods have been employed at various occasions, such as for one-phase mod- 
els in the presence of Griffiths singularities or for the random field Ising model 
[9, 10, 12, 18, 21, 28]. In contrast to the usual case of cluster expansions one does 
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not obtain analyticity (which may not even be valid). In our approach, we loosely 
follow the ideas of Frohlich and Imbrie [21]. For other recent work developing 
their ideas, see [2, 3]. 

5. Cluster expansion of balanced contours 

In this section we perform the zeroth step of the multi-scale analysis for the 
polymer model (19), and set up the cluster expansion for a class of contours the 
weight of which is sufficiently damped. As a result, an interacting polymer model 
is obtained that will be dealt with in the next section. 

Let an integer Iq be fixed. It is supposed to be large enough and the precise 
conditions will be specified throughout the sequel. It plays the role of an rj- 
independent 'cut-off scale'. Given any boundary condition t] (fixed throughout 
this section), we start by defining the set of contours that allow for the cluster 
expansion. Obviously, every bulk contour T has the weight p^{T) = exp(— 2/3|r|). 
For boundary contours, there is no such exponential bound with a strictly positive 
rate, uniformly in rj. Instead, we segregate an 77-dependent subset of sufficiently 
damped boundary contours as follows. 

Definition 5.1. Given rj e O., a boundary contour F is called balanced (or rj-balanced) 
whenever 

E v.>-{i-r)n (23) 

I'd 

Otherwise T is called unbalanced. 

Aset B c 9A is called unbalanced if there exists an unbalanced contour V,d~V — B. 

While the case of large boundary contours will be discussed separately in the 
next section, some basic properties of unbalanced small boundary contours are 
collected in the following lemma. We define the height of any simple boundary 
contour F as 

h{T) = m.&^d[y\dT] (24) 

In order to extend this definition to small boundary contours F such that dV con- 
tains an (exactly one) corner, we make the following construction. If dV is a con- 
nected subset of the boundary with the endpoints [±(A + 1/2), a] and [6, ±(A + 
1/2)], then we define the set -R(F) c {I?)* as the (unique) rectangle such that 
[±{N + 1/2), a], [h, ±{N + 1/2)], and [±{N + 1/2), ±{N + 1/2)] are three of its cor- 
ners. Now the height is the maximal distance of a point in the contour to this 
rectangle, 

/i(F) =maxd[y*,i?(F)] (25) 
The situation is illustrated in Figure 6. 
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Lemma 5.2. Let F be an unbalanced small boundary contour. Then, 

i) E.ear^. <-(!-§) I^ri. 

ii) \dr\ >loh{r). In particular, ifV is simple then \dr\ >lo. 

Proof. For any unbalanced contour F, Definition 5.1 together with the bound |F| > 
\dr\ valid for any small boundary contour implies the inequalities 

-\d-T\< J2 r;.<-(l-f)|r|<-(l-f)|aF| (26) 

Hence, \d^T\ < j^\dT\ and we obtain 

5] < -(1 - hdr\ + |9+F| < -(1 - h\dr\ (27) 
.ear 

proving i). 

If F is simple, then we use (26) again together with the refined relation |F| > 

\dT\+2h{T) to get 

\dr\ > (1 - h\r\ > (1 - i)(|9F| +2/i(F)) (28) 

which implies |c?F| > /o/i(F) > Iq, assuming lo>2 and using that /i(F) > 1 for any 
simple small boundary contour. 

Since the definition of the height is such that the inequality |F| > \dr\ + 2/i(F) 
remains true as well for any small boundary contour F such that dV contains a 
comer, the lemma is proven. □ 
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The union of the set of all bulk contours and of the set of all balanced boundary 
contours is denoted by /Cg. We also write the set of all compatible families 

of contours from /Cq, and for the set of all compatible families of contours from 
ICa \ ATq. Later we will show that, for almost every rj, all large boundary contours 
(i.e. those containing at least one interface) are balanced for all but finitely many 
squares A(A^). 

Formally, the partition function (19) can be partially computed by summing 
over all contours from the set ATq. We start by rewriting partition function (19) as 

Zl^e^p{-E^m J2 Up'i^) H n P'i^') (29) 

Here, the first sum runs over all compatible families d of contours not belonging 
to /Cq, while the second one is over all collections of contours from /Cq, compatible 
with d. Let €^ denote the set of all clusters of contours from /Cq. Then, the cluster 
expansion reads, see Appendix A, 

Zl = exp (-^"(0)) J2 n ^'(r) exp ( E ^oiC)) (30) 

deS^lo red cee^, 

where the sum runs over all clusters of contours from /Cq that are compatible with 
d, and we have denoted the weight of a cluster C by 02(C). Note that the cluster 
expansion was applied only formally here and it needs to be justified by providing 
bounds on the cluster weights. This is done in Proposition 5.4 below. 

Hence, we rewrite the model with the partition function as an effective 
model upon the contour ensemble /Ca \ JCq, with a contour interaction mediated 
by the clusters: 

Zl ^Zl exp (-£;^(0) + 4>1{C)) (31) 

where 

de&la cecg Ted 

After establishing an exponential upper bound on the number of incompatible 
contours in the next lemma, a bound on the cluster weights immediately follows 
by recalling the basic result on the convergence of the cluster expansions [29]. 



KING WITH RANDOM BOUNDARY CONDITION 17 

Lemma 5.3. There exists a constant ci > (independent o/Iq) such that the number of 
all contours V e /Ca, |r'| = n, F' 7^ F is upper-bounded by \r\ e'^^'^,for any F e /Ca and 
n = l,2,... 

Proof. Note that F is not necessarily a connected set. However, the relation F' -t^- F 
implies (F' U ^F') n (F U ^F) ^ 0, and using that F U ^F is connected, we get: 

#{F' : F' 90 F, |F'| = n} < |F U dT\ sup #{F' : e F' U dV , \T'\ = n} 

X* 

< 3|F| sup #{A c (Z^)* connected, x* e A, \A\ < 3n} 

< |F| • 4^"+-^ ^ |r| e*^^" 

by choosing ci large enough. □ 

Assigning to any cluster C e the domain Dom(C) = dC_ where dC — Urec^r 
is the boundary of C, and the length \ C\ — Y^y&c \^\> have the following result. 

Proposition 5.4. There are constants ^4,02 > (independent of Iq) such that for any 
P > ^0/34/ one has the upper bound 

.2/3 



sup 

X* 



E I'^olexp (/-C2)|C| <1 (33) 



x*ec 



uniformly in A. 

Moreover, 4>q{C) only depends on the restriction oft] to the set Dom(C). 

Proof Using Definition 5.1 and equation (21), we have p''(F) < exp(— |^|F|) for 
any balanced contour F. In combination with Lemma 5.3, we get 

5^ \p\T)\ exp[(/ - C2 + 1)|F||] < 5^exp[-(c2 - c, - l)n] < 1 (34) 

reK^ n=l 

provided that C2 is chosen large enough. The proposition now follows by applying 
Proposition A.2, with (^4 = f. □ 

6. Absence of large boundary contours 

By the construction, all unbalanced contours are boundary contours, either small 
or large. In this section we show that unbalanced large boundary contours actu- 
ally do not exist under a typical realization of the boundary condition. This ob- 
servation will allow us to restrict our multi-scale analysis entirely to the class of 
small boundary contours. 
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Lemma 6.1. There is a constant C3 > such that for any N e N and any unbalanced 
large boundary contour F e /CA(iv), the inequality 

5^ r?, < -czN (35) 
x&dr 

holds true. 

Proof. Using the geometrical inequality |r| > 2N + l^+Fj and Dejfinition (5.1), we 
have 



i:,.<-(i-l)|r|.E 

xedT " xed+V 



Vx 



<_(l_l)(2Ar+|a+r|) + |5+r| (36) 
< -2iV(l - h 

where in the last inequality we used that |i9+r| < \dT\ < AN. □ 

Proposition 6.2. There is a constant C4 > such that for any iV G N, 

P{3r e /CA(Ar) large unbalanced} < exp(— C4A^) (37) 

Proof li B (Z dA{N) is a connected set containing exactly two corners, then, using 
Lemma 6.1, 



P{3r e /Ca(jv) large unbalanced : OF ^ B} < P{^r]^ < -C3N} 



„2 



<P{J2vx<-'^\B\}<eM-"fB\) 

xeB 



(38) 



Hence, 

P{3r e /CA(Ar) large unbalanced} 

< P{3r e )Ca{n) large unbalanced : OF ^ B} 

BcdA (39) 

< 2J 8^ exp( — ^/) < — ^exp( — ^) < exp(-C4A^) 

1>2N 

by choosing C4 large enough. □ 

Corollary 6.3. There exists a set Vt* c Vt, P{Vl*} = 1 and a function N* : Q* N 
such that for any b.c. r] e fl* and any volume A = A(A^), > N*{rj), all large boundary 
contours are balanced. 
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Proof. Since 

P{3r e /Ca(jv) large unbalanced} < oo (40) 

N 

the Borel-Cantelli lemma implies 

P{yNo e N : 3A^ > A^o : 3r e }Ca{n) large unbalanced} = (41) 

proving the statement. □ 

We are now ready to prove the almost sure absence of interfaces in the large- 
volume limit. 

Proof of Theorem 4.1. Let rj e Q* n (-0.*) and N > N*{r]). Then, any large bound- 
ary contour T is both r]- and (— ?7)-balanced and, using the Peierls inequality and 
(15), the Gibbs probability of any collection of (possibly large boundary) contours 
Fi, . . . , Fjn, m = 1,2, . . . has the upper bound 

/^A(iv)(ri, • • • , r,n) < ^ max^^ i^I^jv)(ri, • • • , Tm) 

^ 2B (42) 

^ np"'(rO<exp(-/5:ir.|) 

1=1 1=1 

Hence, using Lemma 5.3 and the bound |r| >2N for any large boundary contour 
r, we get 

oo ^ 

^^A{N)i^ ^ large boundary contour) < ^ — - ^ I^a(n) 

(ri,..., 

m=l ' r^,...,r,„ large 
Vi:rjn9A^0 



^E^, E E exp(-H^Eir.i) 

m=l Xl,...,Xm&dAri3xi,...,rm3Xm «=1 (43) 

<-p(~^)Eij(« E ^"^")"' 

° m=l ■ rsa: 
|r|>2]v 

to 

provided that /3 is large enough. Since P{Q* n (— fi*)} = 1, the theorem is proven. 

□ 

As a consequence, all interfaces get P-a.s. and for all but finitely many vol- 
umes uniformly exponentially damped weights. Hence, their Gibbs probabili- 
ties become exponentially small as functions of the size of the system and, there- 
fore, no interfacial infinite-volume Gibbs measure occurs as a limit point, with 
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P-probability 1. While such a result is not sensational in d = 2 (in this case, no 
translationally non-invariant Gibbs measure exists by [1, 26]), similar arguments 
are expected to apply in higher dimensions. 

In the next sections, a perturbation technique is developed that allows us to 
address the question whether non-trivial mixtures of //+ and can occur as limit 
measures. 

7. Classification of unbalanced contours 

We now consider the interacting contour model introduced by the partition 
function (32), defined on the set of unbalanced contours /Ca\A^o- ^ consequence 
of Corollary 6.3, we can restrict our analysis to the set il* of boundary conditions 
under which the set /Ca \ /Cq of unbalanced contours contains only small boundary 
contours, both simple and corner ones. 

Our multi-scale analysis consists of a sequential expansion of groups of unbal- 
anced contours that are far enough from each other. The groups are supposed to 
be typically sufficiently rarely distributed, so that the partition function (31) can 
be expanded around the product over the partition functions computed within 
these groups only. Under the condition that the density of the groups decays fast 
enough with their space extension, one can arrive at an expansion that essentially 
shares the locality features of the usual cluster expansion, at least for P-typical 
boundary conditions rj. To make this strategy work, we define a suitable decom- 
position of the set /Ca\/Cq into disjoint groups associated with a hierarchy of length 
scales. Also, the unbalanced contours close enough to any of the four comers will 
be dealt with differently and expanded in the end. 

Definition 7.1. Assuming Iq to be fixed, we define the two sequences {ln)n=\,2,... und 
{Ln)n=i,2,... by the following recurrence relations: 

i^n = ^, ^n = exp(^) n = l,2,... (44) 

For any n = 1,2,..., any pair of contours F, F' is called L„-connected, if d\r, F'] < 
Ln. Furthermore, fixing a positive constant e > 0, we introduce the A/^-dependent 
length scale 

Zoo = (log TV) (45) 

Introducing the boundary for any set of contours A c /Ca by (9A = \Jre/^dV, 
we consider the T^-dependent decomposition of the set of contours /Ca \ /Cq defined 
by induction as follows. 

Definition 7.2. 1) A maximal Li-connected subset A c /Ca \ /Cq is called a 1- 
aggregate whenever i) \dA\con < k, H) there is no corner such that 
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in.aXy*^gAd[y*,XQ^] < Zoo- use the notation (/C^ q,) /or the collection of all 
'^-aggregates, and write ICI — Uq/Ci q,. 

n) Assume the sets (/CjQ,)j=i,...,n-i have been defined. Then, the n-aggregates are de- 
fined as maximal Ln-connected subsets A c /Ca \ Uj<„/Cj satisfying i) \dA\con < 
/„, ii) there is no corner x^.^ such that max.y*£aA d[y*, Xq^] < loo- The set of all 
n-aggregates is denoted by (/C^ and ICl = UalC'Ji^a- 
To each n-aggregate IC^^ we assign the domain 

^om^lClJ := KgM; d[x\dlClJ<L^} (46) 

Obviously, the set /C^ := /Ca \ {KJl \J KJ^D . . .) need not be empty, and since all 
large boundary contours are balanced, for every contour V e /C^ there is exactly 
one corner such that m.so^y*^aYd[y*,XQ^^ < loo- Hence, there is a natural de- 
composition of the set /C^ into at most four corner aggregates, /C^ = Uj/C^ j, each 
of them consisting of contours within the logarithmic neighborhood of one of the 
corners. In general, any corner aggregate contains both simple and corner bound- 
ary contours. Later we will show that with P-probability 1, every unbalanced 
comer boundary contour belongs to a comer aggregate. In other words, every 
n-aggregate, n = 1,2,... contains only simple boundary contours. 

Remark 7.3. By Definition 7.2, any n-aggregate has a distance at least Ln from all m- 
aggregates, m > n. In this way, in the n-th step of our expansion, after having re- 
moved all lower-order aggregates, we will be able to use the 'essential inde:pendence' of all 
n-aggregates. Namely, on the assumption that Ln is big enough, depending on the ag- 
gregate size In, both the interaction among the n-aggregates and the interaction between 
n-aggregates and m-aggregates, m>n will be controlled by a cluster expansion. 

Our first observation is a local property of the above construction, which will be 
cmcial to keep the dependence of expansion terms to be defined later depending 
only on a sufficiently small set of boundary spins. 

Lemma 7.4. Let a set of small boundary contours A be fixed and assume that 77, r/' e Q 
are such that 77Dom(A) = ^Dom(A)- Then, A is an n-aggregate w.r.t. the boundary condition 
77 if and only if it is an n-aggregate w.r.t. rj'. 

The super-exponential growth of the scales Z„ will imply an exponential decay 
of the probability for an n-aggregate to occur. An upper bound on this probability 
is stated in the following proposition, the proof of which is given in Section 11.1. 

Proposition 7.5. There is a constant C5 > (independent of Iq) such that for any n = 
1,2,... and any connected set B c 9A, 

P{^K:i,a ■■ Con{dlClJ = B}< e-^l^l (47) 
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uniformly in A. 

Note that, given a connected set B c dh., there is at most one aggregate /C^ 
n = 1, 2, . . . such that Con(9/C^ „) = B. 

Corollary 7.6. T/zere exists n** c i^*, P{fi**} = 1 and N** : cu** ^ N, N**{cu) > 
N*{u) such that for any uj e Q,** and any A = A(A'"), N > N**{uj) every aggregate ICl^^, 
n — 1, 2, . . . satisfies the inequality {dJC^ Jcon < ^oo- In particular: 

i) The set Con(/C^^) is a boundary interval and there is at most one corner x}. ^ such 

thatd[x*c^,^,dicij<l^. 

ii) All contours F G /C^ ^ are simple boundary contours. 

Proof Using Proposition 7.5, the probability for any aggregate to occur can be 

estimated as 

Pi^J^la, n = l,2,...: |a/C;^,Jcon > /oo} < Yl ^i^>^n,- ^ Con{IClJ = B} 

B(ZdA conn. 
lB|>(logiV)l + « 

< \dA\ ^ }^N^-o.{lo^Nr ^ ^^jy-S^ 

i>(iogAr)i+« 

(48) 

for any (arbitrarily large) 6 > 0. Hence, 

oo 

J2 P{^K,a, n = 1, 2, . . . : |9/C^,Jcon > /oo} < OO (49) 

N=l 

and the statement follows by a Borel-Cantelli argument. □ 

For convenience, let us summarize the results of the last three sections by re- 
viewing all types of contours again together with their balancedness properties. 
For any r] e fl** and A = A(A^), N > N**{u!), any configuration of contours d e 
possibly contains 

i) Bulk contours (trivially balanced). 

ii) Large boundary contours that are balanced. 

iii) Corner boundary contours that are either balanced or elements of corner ag- 
gregates. 

iv) Simple boundary contours which are balanced or elements of either n-aggre- 
gates, n— 1, 2, . . ., or of comer aggregates. 

8. Sequential expansion of unbalanced contours 

The next step in our strategy is to proceed by induction in the order of aggre- 
gates, rewriting at each step the interacting polymer model (32) as an effective 



KING WITH RANDOM BOUNDARY CONDITION 



23 



model over the contour ensembles /Ca \ (/Cg U /C?), /Ca \ (/Cg U /C? U K,^), etc. At 
the n-th step, a compatible set of contours inside all corner and all normal m- 
aggregates, m > n, is fixed, and we perform the summation over contours in all 
normal n-aggregates. This is a constrained partition function which is approxi- 
mately a product over the normal n-aggregates. By the construction, the latter 
are sufficiently isolated on the scale L„, which will allow for the control of the 
remaining interaction by means of a cluster expansion. At the end, we arrive at 
an effective model over the contour ensemble /C^, which is the union of (at most 
four) corner aggregates. In large volumes, the corner aggregates become essen- 
tially independent, the error being exponentially small in the size of the volume. 
The reason we distinguish between the n-aggregates and the comer aggregates 
is that the partition function within the former allows for a much better control, 
which will be essential in our analysis of the characteristic function of the random 
free energy difference log — log Z^^ in Section 10. Note that the lack of detailed 
control around the corners is to be expected as there may more easily occur some 
low-energy (unbalanced) boundary contours, but at most of logarithmic size in N. 
The n-th step of the expansion, n > 1, starts from the partition function. 



which in the case n = 1 coincides with (32). Here, ^>„_i is the set of all compatible 
families of contours from /C^„_i := /C](\ (/CqUAI^'U. . .UlCl_i), i.e. with all normal m- 
aggregates, m < n — 1, being removed. Furthermore, we use the notation ^^-i 
the set of all (n— l)-clusters. Here, the 0-clusters have been introduced in Section 5, 
and the clusters of higher order will be defined inductively in the sequel. 

In order to analyze partition function (50), we follow the ideas of Frohlich and 
Imbrie [21], however, we choose to present them in a slightly different way. Ob- 
serving that, by construction, the family of aggregates compose a 'sparse set', one 
is tempted to approximate the partition function by a product over the aggregates 
and to control the error by means of a cluster expansion. However, to make this 
strategy work, we need to 'renormalize' suitably the contour weights. Namely, 
only the clusters that intersect at least two distinct aggregates generate an inter- 
action between them, and are sufficiently damped by using the sparsity of the 
set of aggregates. On the other hand, the (sufficiently short) clusters intersecting 
a single aggregate cannot be expanded, and they modify the weights of contour 
configurations within the aggregate. An important feature of this procedure is 
that the weight of these contour configurations is kept positive. In some sense, it 
is this very renormalization of the weights within each aggregate that can hardly 




(50) 
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be done via a single expansion and requires an inductive approach. In what fol- 
lows, we present this strategy in detail, via a number of steps. 

8.1. Renormalization of contour weights. For any compatible set of contours 9 C 
/C^, define the renormalized weight 

pna) = exp(- Yl '^n-i(^^))n^'(n (51) 



ceel_-^ red 



Note that the above sum only includes the clusters of length smaller than L„. 
By construction, any such cluster is incompatible with at most one n-aggregate. 
Hence, the renormalized weight p''(r) factorizes over the n-aggregates and we 
have p^ip) = riaP^'l^") where 9" = 9 n /C^^^- Therefore, formula (50) gets the 
form 



ae^^„ red a^e^^ 

(C*a)V(Co<.a»; \C\>Ln) 

J2 n'«'(nexp(- Y '^n-l(C)) 



de&>„ red ceel_-^ 



(52) 



X 



J2 P'{dnew{- E ^n-iiC)) 



d"e^Z ceel_-^;\c\>Lr, 

Defining the renormalized partition function of the contour ensemble /C^'" 

as 

4v= E '^'(^") (53) 

d^e&loc 

and using the shorthand 

9") = ct>:-iiC)l{c.9n; \c\>L„} (54) 

we obtain 



desil^ red cec 



V 
n-l 



(55) 



7^ 
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where 9"'" = 9" n /C^ is the restriction of 9" to the n-aggregate /C^ q,- I^i the 
last expression, the second sum contains the interaction between n-aggregates, to 
make a correction to the product over the renormalized partition functions Z^^. 



8.2. Cluster expansion of the interaction between n-aggregates. Now we em- 
ploy a trick familiar from the theory of high-temperature (Mayer) expansions, and 
assign to any family C C of (n — l)-clusters the weight 

<iC) = P'i^n n (e-^"-^^'^") - 1) (56) 

lia^A dn^^v CeC 

See Figure 7 for an example of a family of 1 -clusters that generically gets a non- 
trivial weight according to this construction. 

Definition 8.1. Any pair of (n — l)-clusters Ci, C2 G is called n-incompatible, 

n 

Ci C2, whenever there exists an n-aggregate IC^^^ such that Ci ^ IC^^^ and C2 ^n,a- 
In general, the sets Ci,C2 C are n-incompatible if there are Ci e Ci, C2 e C2, 

n 

Ci ^ C2. 

One easily checks the following properties of the weight tyJJ(C). 
Lemma 8.2. For any set of {n — l)-clusters C e 

i) supJ<(C)|<nc.c(el'^"-^''^l-l). 

ii) IfC = Ci U C2 such that Ci A C2, then w'^iC) = w:^{Ci) ^(Cs). 

iii) The weight w^KC) depends only on the restriction ofrj to the set (Ucec Dom(C)) U 
(U^ Dom(/C^ ^) where the second union is over all n-aggregates IC^ ^ such that 

In the second sum in (55) we recognize the partition function of a polymer 
model with the polymers being defined as the n-connected subsets of which 
are incompatible if and only if they are n-incompatible. Treating this polymer 
model by the cluster expansion, and using the symbols for the set of all clus- 
ters in this polymer model and i/j^.iD) for the weight of a cluster D e D'!l , we 
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Figure 7. A pair of 2-compatible families of 1-clusters Ci,C2 C 
intersecting 2-aggregates K^^^, a = 1,2, 3, 4. The dashed rectangles 
illustrate 1-aggregates which have become parts of the 1-clusters. By 
construction, w^(Ci U C2) = w^(Ci) wl{C2). 



get 



= exp( E E n^^(r) (57) 

xexp(- E ^l-AC)- E 



Defining the set of all n-clusters <L1 = C^-i U and the weight of any n-cluster 

c e c;^ as 

we finish the inductive step by obtaining the final expression 

Zl = n 4V exp ( E rn{D)) (59) 
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with the partition function of a new interacting polymer model 

^n\i = E ^^p(" ^ '^^^^O n/''(r) (60) 



de^l„ ceel red 

^ c*a 



We need to extend the notion of domain from the set of (n — l)-clusters ^^-i to 
the set of n-clusters C^. Realizing that any n-cluster D e D^is a collection (Cj) of 
L„-connected families of (n — l)-clusters, Ci = {C-), we first introduce the domain 
of any such family Ci as Dom(Ci) = Dom(C/). Next, we define 

Dom(D) - U Dom(C,) U (J Dom(/C:i,J (61) 
Furthermore, the length of the cluster is defined as 

i^i-Eic^i = EEi^^i (62) 

i is 

Note that this is possibly much smaller than the diameter of the cluster, since the 
sizes of the n-aggregates in the domain of D are not counted in the length of the 
cluster. The reason for this definition is that the cluster weights are not expected 
to be exponentially damped with the cluster diameter. Note, however, that the 
probability of a cluster to occur is exponentially damped with the size of the n- 
aggregates in its domain. 

In the next proposition, we provide uniform bounds on the n-cluster weights. 
For the proof, see Section 11. 

Proposition 8.3. There is jS^ > such that for any (5 > Iq(5^, rj e fl**, A = A(A^), 
N > N** (rj), the inequalities 

sup E exp(f|i?|)|V'OTI<2-" n=l,2,... (63) 



Devi 



and 



supsup E cxp(f \ct>l{C)\ < 1 (64) 

'0 



X* 

x*ec 



hold true. 

Moreover, ifC e and 77'|Dom(c) = ^lDom(c)/ ih^n also C e €^ and (p''J^{C) = (pliC). 
Similarly, D e Vy^andr]'\Bom{D) = v\Bom{D) implies both D e 'D'^ and il}'!^ {D) = ipli{D). 
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8.3. Expansion of comer aggregates. For any finite square A = A(A'^) and rj e Q, 

all aggregates from the set U„/C'; are expanded in a finite number of steps. After- 
wards, all comer aggregates are treated by a similar procedure. Throughout this 
section, we use the notation uq for the highest order in the collection of all normal 
aggregates. The expansion goes similarly as in the case of normal aggregates, so 
we only sketch it. 

The renormalized weight of any compatible family of contours d C /C^ is de- 
fined by the formula 

pn9)=exp(- Yl <^no-i(C))n'«'(n (65) 

which factorizes over the comers, p''((9) = Yl- p^{d fi /C^ J, assuming A(A^) to be 
large enough. Clusters Ci, C2 C are called 00-incompatible whenever there is 
a corner aggregate /C^ ^ such that Ci /Coo,? and C2 /Coo,i. Defining the weight 

w\C) = — i^p Yl P"^^) n (^^-^-^"^'"^ - 1) (66) 

where 
and 

0^(C,9) = C(C)l{c^a;|ci>2;.} (68) 

an obvious variant of Lemma 8.2 holds true and w'^{C) factorizes into a product 
over maximal connected components of C w.r.t. 00-incompatibility. Treating these 

as polymers in a new polymer model with cxD-incompatibility used as the incom- 
patibility relation, and using the notation for the set of all clusters in this 
polymer model and 'ijj2^{D) for the cluster weights, we obtain as the final step of 
the sequential expansion, 

^no+l - e^p( E ^-(^)) n ^ki (69) 

Proposition 8.4. There exist constants (3q > /35,cq > such that for any (3 > IqI3q, 
77 e Q**, and volume A{N), N > N**(r]), one has the bound 

E \ro.m<e-^'- (70) 
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Gathering all expansion steps, we arrive at the final expression for the partition 



function Z\ in the form 



cee^ n>i DeD'^ De^2o 



(71) 



n>l a 



The terms collected on the first line contain the 'vacuum' energy under the bound- 
ary condition f], together with the contributions of clusters of all orders. Recall 
that the latter allow for a uniform exponential upper bound. On the second line 
there are the partition functions of all n- and all corner aggregates. Although we 
can provide only rough upper bounds for these terms, a crucial property to be 
used is that the probability of an aggregate to occur is exponentially small in the 
size of its boundary, see Section 7. In this sense, the above expansion is a natural 
generalization of the familiar 'uniform' cluster expansion [29]. 

8.4. Estimates on the aggregate partition functions. In expression (71) we do not 
attempt to perform any detailed expansion of the aggregate's (log-)partition func- 
tions and • via a series of local and exponentially damped terms. Instead, 
we follow the idea that a locally ill-behaving boundary condition forces a partial 
coarse-graining represented above via the framework of aggregates of different 
orders. Although the detailed (cluster expansion-type) control within the aggre- 
gates is lost, we still can provide generic upper bounds on these partition func- 
tions. Notice a basic difference between n-aggregates and comer aggregates: The 
former contain only simple boundary contours the weights of which exponen- 
tially decay with the height of the contours. In some sense, the partition functions 
Z^^a c^ri be compared with the partition function of a Id interface to get an up- 
per bound. On the other hand, the corner aggregates are ensembles of contours 
the weight of which obey no uniform exponential bound with the space exten- 
sion of the contours, and allow possibly for a non-trivial 'degeneracy of vacuum'. 
As a consequence, only rough (counting-type) estimates can be provided for the 
partition functions Z^^^. 

Lemma 8.5. There are constants 07,07 > (cr I if P ^ 00) such that for any n- 
aggregate IC^^^, one has the hound 

\ogZl,<c,\dKZJ (72) 

For any corner aggregate KP^ ^, 

log Zl,<c',ll (73) 
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9. Asymptotic triviality of the constrained Gibbs measure v\ 

As the first application of expansion (71) we prove that the weak limit of the 
constrained measure v\ coincides with the '+' phase Gibbs measure finishing 
the first part of our program. 

Proposition 9.1. There exists a constant c > such that for any (5 > Iq(5q (with the (5^ 
the same as in Proposition 8.4), any r) e fi**, and X c 7? finite, 

II^Iw-/^lU = 0(e-"^) (74) 
In particular, liniTv^oo '^A(Ar) = A*^' P-a.s. 

Proof. The idea of the proof is to express the expectation i^^^^if) of any local 
function / as the sum of a convergent series by using the multi-scale scheme 
developed in the last section, and to compare the series with a standard cluster 
expansion for i^]^^y The difference between both series is given in terms of clus- 
ters both touching the boundary and the dependence set of /. Restricting only to 
the boundary conditions rj e fl** and volumes ^{N), N > N**{i]) and using the 
exponential decay of the cluster weights, we prove the exponential convergence 

For notational simplicity, we only restrict to a special case and give a proof of 
the equality 

lim uliao = -1) = /x+((7o = -1) (75) 

The general case goes along the same lines. 

Assuming a e Qj, observe that ao — —Hi and only if the set T>A{a) contains an 
odd number of contours F such that G Int(F). In an analogy with (19), we write 
the ^^-probability that o-q = — 1 in the form 

<i^o = -1) = 4 E ^a(\^) n ^'(r) (76) 

A AcA reA 

where we have used the shorthand A iz A for any compatible family of contours 
in A such that card (A) is an odd integer and e Int F for every F e A . Further- 
more, Z]((\A) is the partition function 

des>A{\A) red 

of a polymer model over the restricted ensemble /Ca(\A) c /Ca of all contours 
F such that i) F ~ A, and ii) ^ Int(F). We can now repeat the same proce- 
dure as in the last sections, but with the contour ensemble /Ca being replaced by 
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X^''(\A). A crucial observation is that all contours from the set /C \ /^(yA) are bal- 
anced, at least for all rj G Vl** and provided that the volume A(A^) is large enough. 
Hence, the sets of unbalanced contours coincide for both contour ensembles /C 
and /C''(\A), hence, the same is true for the collections of both n- and comer ag- 
gregates. Finally, we compare the terms in the expansions for Z\ and Z^(\A), and 
arrive at the formula 

^ Ce<t^\<t^{\A) n>l D&^l\^l{\A) DGS)2o\S2o(\A) 

(78) 

where each of the three sums runs over all (0-, n-, or oo-)clusters that are either 
incompatible with A or contain a contour F, e Int(r). By construction, each 
n-, respectively oo-cluster is further required to be incompatible with an n-, re- 
spectively corner aggregate, and since their weights are uniformly exponentially 
bounded by Propositions 8.3-8.4, we get the uniform upper bound 

2im 



sup 

A 



log- 



7"^ 



<c\A\ (79) 



with a constant c large enough, as well as the existence of the limit 

limlog5^ = -^0o(^^) (80) 
^ c 

where the sum runs over all finite 0-clusters in that are either incompatible with 
A or contain a contour surrounding the origin. 

Since every F G A surrounds the origin, it is necessarily balanced and satisfies 
p'?(r) < exp(-|^|r|). Combined with (79)-(80), one easily checks that 

limul{ao = -l)= Yc^{-Y,^PoiC))l[piT) (81) 
Acz2 c reA 

and the convergence is exponentially fast. Obviously, the right-hand side coin- 
cides with the limit limjs_ iJ,^~^^ {ao — —1) = //■^(ctq = —1), which finishes the 
proof. □ 

10. Random free energy difference 

In this section we analyze the limit behavior of the sequence of the random free 
energy differences 

F]l = log Zl- log ZJ^"^ (82) 
In order to show that the probability that takes a value in a fixed finite inter- 
val is bounded as 0(A^~2+") with a > 0, we can use the local central limit upper 
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bound proven in Appendix B, provided that a Gaussian-type upper bound on the 
characteristic functions of the random variables can be estabHshed. The basic 
idea is to prove the latter by employing the sequential expansion for log Z''l devel- 
oped in section 8 and by computing the characteristic functions in a neighborhood 
of the origin via a Mayer expansion. However, a technical problem arises here due 
to the high probability of the presence of comer aggregates. That is why we need 
to split our procedure in two steps that can be described as follows. 

In the first step, we fix the boundary condition in the logarithmic neighborhood 
of the corners and consider the random free energy difference conditioned on 
the fixed configurations. For this conditioned quantity a Gaussian upper bound 
on the characteristic function can be proven, implying a bound on the probability 
that the conditioned free energy difference P-a.s. takes a value in a scaled interval 
{aN^, bN^). This can be combined with a Borel-Cantelli argument to exclude all 
values in such an interval, at least P-a.s. and for all but finitely many volumes 
from a sparse enough sequence of volumes. 

In the second step, we consider the contribution to the free energy difference 
coming from the corner aggregates. However, their contribution to the free energy 
will be argued to be of a smaller order when compared with the contribution of 
the non-comer terms. 

Note that we also include the oo-clusters in the first step. Because we have 
uniform bounds in i] for the oo-cluster weights, we are allowed to do so. 

The free energy difference F^ can be computed by using the sequential expan- 
sion (71). For convenience, we rearrange the terms in the expansion by introduc- 
ing 

U^'iB) = J] J]logi^^^l{Dom(C„)=B} + J]logi^^,l{Dom(K^_J=B} 

n a i 

+ </'2(C)l{Dom(C)=B} + XI XI V'n(^)l{Dom(D)=B} (83) 
C n D 

{Dom(D)=B} 

D 

for any set B c 9A. Note that any function U'^{B) only depends on the restriction 
of 77 to the set B. Using the notation W{B) = U'^{B) - U'^iB), the expansion for 
the free energy difference reads, formally. 



F2 = 2/5 X + E ^'(^) 

x£dA BcdA 



(84) 
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Obviously, no bulk contours contribute to U^{B). Using the notation dAc,i '■— 

{y* e OA : d[y*,X(^^] < 21^} and dAc uf^idAc,i, we consider the decomposi- 
tion Fl^Fl + Fl, where 

and 



5 — Dom(S)(z:aA^ 



xedAc BCdA 
Dom(B)caAc 

The jfirst term, F^{B), can be analyzed by means of the Mayer expansion of its 
characteristic function 

^l{t):=E[exp{ttF]l)\vaAc] = E[exp{2tti3 Yl ^a.) E 11 (^"*''''^''^ " 1) I 

xe dA\dAc B BeB 

B 

(87) 

where we have assigned to any family B of subsets of the boundary the weight 

^ xe dA\dAc BeB (88) 

X l{VBeB:B!Z:aAc} 

and have introduced the notation 

*o(^) = E[exp{2itprio)] = cos 20 (89) 

Observing that 

w{Bi U B2 I r] 9Ac ) = uj{Bi I r] dAc ) w{B2 \ ^ qac ) (90) 

whenever BiH B2 = % for any Bi E Bi and B2 G B2, the last sum in equation (87) 
is a partition function of another polymer model and using the symbols ?B, Q3i, . . . 
for tine clusters in this model and wj for the cluster weights, we get 

nit) = [vl/olt)]!^'^^^^-! exp Y wJiB \ vaAc)] (91) 

B 

A crucial observation is that for any 77 e Q** and A{N), N > N**{r]) no corner ag- 
gregate contributes to the weight wt{B) for any B. On the other hand, the partition 
function of any n-aggregate is balanced by a small probability of the aggregate to 
occur. Another observation is that every weight Wt{B) is of order O(t^) due to the 
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symmetry of the distribution P. To see this explicitly, formula (88) can be cast into 
a more symmetrized form. 



(92) 



Bee 

where Supp(i3) := UegbB and 



VdAc 



j^m.^Psiny if card(S) = 2A; - 1 

■ [cosy if card(,B) = 2A; keN ^ ^ 

In section 11.5 we give a proof of the following upper bound on the corresponding 
cluster weights: 

Lemma 10.1. There exist constants (3s, lo > 0^ such that for any [3 > (Sglo there is 
to = ^o(/5) > 0/or which the following is true. For any tj e Q** and A = A(A^), 
N > N** {rf), the inequality 

sup \^J{B\r^eKa)\<UH'' (94) 

is satisfied for all\t\ < to. 

With the help of the last lemma, it is easy to get an upper bound on ^\{t): 
Lemma 10.2. Under the assumptions of Lemma 10.1, we have 

nit) < exp(^~(3H'\dA{N) \ aAc(iV)l) (95) 
for all \t\ <to,v ^ ^**, and N > N**{r]). 

Proof. It immediately follows by combining Lemma 10.1, equation (91), and the 
bound *o(^) < ew[-(3H% □ 

For the comer part of the free energy difference we use the next immediate 
upper bound: 

Lemma 10.3. Given rj e ft** and (5 > (5qIq, then F^^^^^ = 0{N^)for any 5 > 0. 

Proof. Using Proposition 8.4 and Lemma 8.5, we have YIib^okc ^ ^iPoo) 

and the above claim immediately follows. □ 



^Recall that the construction of aggregates depends on the choice of Iq. 
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Proof of Proposition 3.5. Combining Lemma 10.2 with Proposition B.l in the appen- 
dix, we get 




for any a,T > 0. By Lemma 10.3, F can be replaced with the full free energy 
difference F. As a consequence, 

^ 7V5-P{|f;^ |<r}<oo (97) 
and the proof is finished by applying Proposition 9.1. □ 

11. Proofs 

In this section, we collect the proofs omitted throughout the main text. 

11.1. Proof of Proposition 7.5. In order to get the claimed exponential upper 
bound on the probability for an n-aggregate to occur, we need to analyze the way 
how the aggregates are constructed in more detail. We start with an extension of 
Definition 7.2. Throughout the section, a finite volume A = A(A'^) is supposed to 
be fixed. 

Definition 11.1. For every n = 1,2,..., any maximal Ln-connected subset A c /C \ 
(/Cq U /Ci U . . . U /C^_i) is called an n-pre-aggregate. 

Obviously, n-aggregates are exactly those n-pre-aggregates A that satisfy the 
condition |9A|con < In- Moreover, every n-pre-aggregate can equivalently be con- 
structed inductively by gluing pre-aggregates of lower orders: 

Lemma 11.2. Every n-pre-aggregate A„ is the union of a family of {n—l)-pre-aggregates, 

An = UaA°„i. Moreover, 

i) Each {n — l)-pre-aggregate satisfies |9A^_i|co« > In-i, 

ii) The family (A"_ Jq, is Ln-connected. 

Proof. For n = 1 the statement is trivial. 

Assume that n > 2, and let A„ be an n-pre-aggregate and F e A„ be a contour. 
Then, there exists an {n — l)-pre-aggregate A"_i such that F e A"_i (otherwise F 
would be an element of a fc-aggregate, k < n — 2). Moreover, since A"_j^ is not an 
(n — l)-aggregate by assumption, it satisfies |(?A^_^|con > ^n-i/ proving i). 
The claim ii) is obvious. □ 

Lemma 11.3. Let A hy any family of unbalanced contours. Then, 

i) There exists a subset A c A such that 
a) dK = dA, 
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b) ifTi, Fs e A are any three mutually different contours, then dVi n dV2 n 
dV^ = 0. 
ii) The inequality 




holds true. 

Proof i) Assume that Fi, F2, F3 c A is a triple of mutually different contours such 
that dTi n 9F2 n 9F3 7^ 0. Since dTi, i = 1,2,3 are cormected subsets of OA, it is 
easy to realize that, up to a possible permutation of the index set {1, 2, 3}, one has 
dTi c aFsUaFg. Hence, 5(A\{Fi}) = dA. Since the set A is finite, a subset A c A 
with the claimed property is constructed by iterating the argument. 

ii) Let A c A be the same as in i). Then, using Lemma 5.2, the inclusion- 
exclusion principle implies 

xedA rezA xedr (r r')r:A xe drndT' 

<-(i-r)El^r|+ V laFn^F'l (99) 

reA (r,r')cA 
<-(l-i)|aA| 

□ 

It remains to prove that one still gets a large deviation upper bound by replacing 
the sum over the boundary sites x e dA in equation (98) with the sum over all 
X G Con(9A), provided that A is a pre-aggregate. Technically, we need to exploit 
the basic feature of any pre-aggregate A that the set Con (9 A) \ dA is not 'too 
big'. A minor complication lies in the fact that the boundary distance d[d^, d^] 
is allowed to exceed the contour distance ^[7, 7']. To overcome this difficulty, it is 
useful to define 

C^{dA) = dAu{xe Con(aA); V7 G A : d[x, ^7] > (100) 

lo 

for which the first equation in the proof of Lemma 5.2 implies the upper bound 

|aA|con<(i + r)lc^(^^)l (101) 

We are now ready to prove the following key estimate from which Proposi- 
tion 7.5 immediately follows by using a large deviation upper bound. 
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Lemma 11.4. Let A be an n-pre-aggregate, n— 1,2,... Then, 

^ -l\dA\,an (102) 

a:e Con(aA) 

uniformly in n. 

Proof. We prove by induction in the order of the pre-aggregates the refined bound 

T ^ 

J2 r/.<-(l-35]-^)|Con(aA)| (103) 

xe Coii(aA) »=i 

for any n-pre-aggregate A, from which the statement follows by using the defini- 
tion 7.1 of length scales Z„ and and equation (101). Indeed, one obtains then 

OO J- 

^ - 3E |c5n(9A)| + l Con(aA) \ c5n(9A)| 

i;eCon(9A) i=l (104) 

1 2jaAU - 1 



+ 1 



< --I9AI, 



First, assume that A is a 1-pre-aggregate, and let A = U^^^Ai be the (unique) 
decomposition of A into disjoint subsets such that dA = U^^dAi is the decom- 
position of dA into maximal connected components. For convenience, we use the 
notation Jj := dAi. Considering furthermore the decomposition Con(9A) \ dA — 
If^IiGk into maximal connected components, the set Con(9A) can be finally writ- 
ten as the union 

Con(9A) = (uriiJ,) U (U^-^G,.) (105) 
of disjoint connected subsets, which satisfy the inequalities \Ji\> Iq and IG^I < Li, 
for alH, /c = 1,2,... Using Lemma 11.3, we have Ylx&Ji Vx < ~ ~ ^) I I' ^^"^ since 



EkJ \Gk\<j^ TZi 1-^.1/ we finally get 



m—l 



V- v^v^ V^i^i^ /i Li +4JCon(5A)| 

x6 Con(aA) »=1 3;eJi fe=l «o ^]^Q^^ 

Q r 

<-(l-T^)|Con(aA)| 
to 

provided that, say, Li > 4. 

Next, we will prove the statement for an arbitrary n-pre-aggregate A. By Lemma 11.2, 
A is the union of a family of (n — 1) -pre-aggregates, A = VJiA\^_^. In order to gener- 
alize our strategy used in the n = 1 case, we consider the (possibly disconnected) 
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boundary sets Jj = Con((9Aj^_ J, and the family of connected sets (G'i)i=i,2,... de- 
fined as the maximal connected components of the set Con(9A) \ Uj Coii(A^_i). 

Note that = #{ J^} - 1 and the identity Con(A) = (U^ J^) U (UjGj). Hence, 

by using the induction hypothesis, 

E •'-EE''. + EiG*isK-3E^) + ^]^|^ 

a;£ Con(aA) «=1 ^eJi fc=l »=1 (n-i ^^^^^ 

^-(l-3EA)|cSn(9A)| 

as required. □ 

11.2. Proof of Proposition 8.3. The proof goes by induction in the order of aggre- 
gates. 

The case n = 1. 

As the initial step we bound the sums over 1-clusters in D'1. Recall that the 1- 
clusters consist of 0-clusters which connect 1 -aggregates /C^ q,. Throughout this 

section we use the shorthand /9 := £. 

From Proposition 5.4 we know that for any integer ro, 

J2 |02(C)|exp(^(2-(l/8))|C|)<l 



-0- 

C3x 



which implies 



J2 103(^)1 exp (2/3(1 - (1/8)) < exp (-/3ro/8) (108) 



CeCg: lCl>ro 
C3x 



We split the procedure into four steps as follows. 

Part 1. For any 1-cluster in D^, none of its 0-clusters contributes to the dressed 
weight of a 1-aggregate. Hence, all these 0-clusters have at least size Li. Moreover, 
they are incompatible with a 1-aggregate /C^ q,. Using Lemma 5.2 and choosing 
ro — Lim (108), this results in the inequality 

J2 |02(C)|exp(2^(l-(l/8))|C|)<Z?exp[-(^Li)/8] <2-2 (109) 
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Part 2. In order to prove the convergence of the cluster expansion resulting from 
the Mayer expansion, we apply Proposition A.2. As our initial estimate, we get, 
using (109) and since 

C^C ^ 3a such that C, C / JCl^ 

the inequality 

J2mC)\ew{2P{l-{l/8))\C\) 

< E E exp(2^(l-(l/8))|C|)|02(C)| (110) 

< ^-Wia ^ Co] 

Part 3. Using Lemma 8.2, the weight of any set of 0-clusters appearing in the 
Mayer expansion is bounded as 

K(C)|<n(el^°^^^'-l)<n2|'^2(C^)l 

Hence, by using Proposition A.2, we obtain the bound 

|V^?(Ci)|exp[(2/3(l - (1/8)) - l/2)\C,\] < 2-^#{/C?,, 7^ Co} (111) 

Taking now Co e €q such that /C^ is the only 1-aggregate satisfying Cq oo JCf, 
inequality (111) yields 

J2 |^?(Ci)|exp[(2^(l - (1/8)) - l/2)|Ci|] < 2"^ (112) 

Part 4. In order to bound the sum over all 1 -clusters Ci e X>i such that Ci 3 x 

and |Ci| > ri, we use that |Ci| > Li and write 

J2 mC,)\e^v [(2^(1 - (1 + l/2)/8) - l/2)|Ci|] 

(113) 

<E E |V'?(Ci)|exp[(2^(l-(l + l/2)/8)-l/2)|Ci|] 
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Substituting (112), we obtain 

(113) < ^2-iexp(-(^/8) •max[d(/C?,,,x),ri]) 

oo 

< XI 2-^exp(-e/g-max[i?,ri]) 
The last sum can be estimated by a partial integration and we finally get 



(113)<exp(-;3ri/8) 



9 16ri 



< 2-^ ■ Arl exp (-/9ri/8) 



where we have used that ri > Li and that Li is large enough. 



Induction step. 

The induction hypothesis reads 



E lV'?(a)|exp 



Ci3x: \Ci\>ri 



i+1 i 

(2;g(i-E(i/2)V8)-E(V2y)ia 



j=0 j=l 

< 4 • 2-V| exp (-^(l/2)^+Vi/8) 



for any 1 < i < n — 1. 

Part 1. As in part 1 of the n—1 case, we want to prove first that 



E l'^n-i(C)|exp 



n-1 



(2^(1 -E(l/2)V8)-E(V2r) 1^1 

j=o j=i 



< 2 



-n-1 
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Recalling Definition (58) for (t>l-i{C), we know that (t)l_^{C) = ^]{C) for any C e 
S)^. Hence, using (115) with r, = L„, we write 



n— 1 



(116) <Z^^ J2 \i^nC^)\ew 



i=l Ci3x: \Ci\>L„ 
71-1 

i=0 



<2-"-i-32/2L2exp -/3(l/2)'^L„/4 



n— 1 



(2^(1 -5](1/2)V8)- 5^(1/2)^) 



n-l 



(-^ J](1/2)V4- J](1/2)^)l„ 



j=i+i 



< 2 



-n— 1 



(117) 



where we have used that Z„ = exp (L„/2"^) and ;5 is large enough. This proves 
inequality (116). 



Part 2. Similarly as in the n — 1 case, we prove by using (116) the inequality 



5] \<f>U{C)\e^p 



C^Co 
n— 1 



n-l 



(2^(1-5](1/2)V8)-5](1/2)^)|C| 



3=0 



< 



^-""-Wla ^ Co} 



Part 3. By construction, any n-cluster C„ e consists of a family of 0-clusters 
Co e ^rid i-clusters Ci e Sg/ < i < n — 1, which are all incompatible with /C^. 
Using Lemma 8.2 again, we have the upper bound 



n—l 



i<(Cn)i<n n 

where we have identified ^o(.) = and Dq = €q. Applying Proposition A.2 

with z{C) = 2\^1'{C)\ then gives 



n n 

(2^(1-5](1/2)V8)-5](1/2)^)|C| 

3=0 j=l 



< 2-"#{/C^,„ 7^ Co} 



Taking again Cq e such that IC^^^ Cq implies the inequality 



E \rn{Cn)\eW 



(2;3(1-E(1/2)V8)-E(1/2)^)|C„| 



3=0 



3=1 



< 2-" 



(118) 
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Part 4. Repeating the argument for the n—1 case, we obtain the inequality 

n+l 



E \rn{Cn)\eW 



Cn3x, \Cn\>rn 



2/3(1- 5^(1/2)78) -5^(l/2y)|C„| 



j=0 j=l 
„2 / /I /o^n+l, 



(119) 



< 2-" • Kexp (-(l/2)"+^/?r„/8) 



Using that r„ > L„ for any e 2)^ and choosing Vn — Ln proves the proposition 
for the weights V'n/ ^ = 1)2, 

Equation (58) reads that (t)l{C) = ip]{C) whenever C e Ti] and j < n. Using 
further that = €oUS)iU---US)^ and summing up the cluster weights of the 
clusters of all orders yields inequality (64), which finishes the proof. 

11.3. Proof of Proposition 8.4. Let uq be the same as in Section 8.3. Due to the 
second part of Proposition 8.3, 

sup Yl \K{C)\exp[{(3/4lo)] \C\ < 2exp (-(3/3/4/o)ro) 

^ C3x, |Cl>ro 

According to the definition of the corner-aggregates, we have 

J2 \<Pl{C)\exp[{P/4lo)\C\] < 24exp(-/oo(3/3/2/o)) < 2"^ exp (-/oo(/5//o)) 

Applying Proposition A.2, we obtain 

E |^^(C)|exp[(/5/4/o)] \C\ < 2-'exp{-U(3/lo)) 

which implies 

|^S^(i^)|<exp(-(3/3/2Zo)/oo) 

11.4. Proof of Lemma 8.5. Let 77 e O.** and /C^^ be an n-aggregate, n = 1,2,.... 
Recall that 
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where 

p\d) = Ylp\T)exp[- Yl ^n-i{C)) (121) 

red ceel_-, 

C'^d; \C\<L„ 

Using the 77-unif orm bounds 

p''(r)<exp[-2/?(|r|-|ar|)] (122) 

and 

sup E l€-i(G)l<exp[-^] (123) 

n— 1 

for all n = 1,2,..., one can subsequently write (for simplicity, we use the short- 
hand £ — exp [—7^] below): 

^n,"< E n^^pH2/5-£)|r| + 2/3|9r|] 

< g(.+4e-^+^)|a;c;^,.| ^ exp[-(2/3 - £)|r| + (2/3 - £ - 4e-2^+^))|9r|] 

< e^^+^^-^'^^^-l E J] J] exp[-(2/3-£)|7| + (2/3-£ -46-2^^+^)1971] 

\dicia\ 



< e(-+4^-^^+^)l^^^-l{l + ^exp[-(2/3 - 5)|7| + (2/3 - £ - 4e-2^+^))|a7|]} 



73p 



(124) 



where the last sum runs over all pre-contours (= connected components of con- 
tours) such that a fixed dual bond p — {x, y)*, d{x, A^) = d{y, A'^) = 1 is an element 
of 7 and it is the leftmost bond with these properties, w.r.t. a fixed orientation on 
the boundary. To estimate this sum, we associate with each pre-contour 7 a path 
(= sequence of bonds; not necessarily unique) starting at p. Every such a path 
consists of steps choosing from three of in total four possible directions. One eas- 
ily realizes that, for every such a path, the total number of steps to the right is 
bounded from below by \d^\. Hence, the last sum in (124) is upper-bounded via 
the summation over all paths started at p, so that to each step going to the right 
(respectively to the left/up/down) one assigns the weight e"^^ ^'^'^^ (respectively 
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g-2/3+£^^ which yields 



5]exp[-(2/3 - £)|7| + (2/3 - £ - 4e-^^+^))|97|] 

(125) 

n=l 

All in all, one obtains 

< e^^+^^-^'^^^l^-l (126) 

proving the first part of the statement. 

The proof of the second part is trivial by counting the number of all configura- 
tions in the square volume with side 21 oo- Note that the latter contains all contours 
r e 9 for any configuration d e ^ and that the weights of all clusters renormal- 
izing the contour weights are sumable due to Proposition 8.3. 



11.5. Proof of Lemma 10.1. Due to Proposition A.2, it is enough to show that the 
inequality 

J2 \MB I VaAc)\ ^Ml^'t'l Supp(B)l) < l^H' (127) 

B:x*eSupp(B) 

holds true for all |i| < to, with a constant to > 0. Remark that the RHS of the last 
equation is not optimal and can be improved, as obvious from the computation 
below. 

In order to prove (127), we use the symmetric representation (92) of the weight 
wt{B I rjdA^), the lower bound ^o(^) > which is true for any a > provided 
that \t\ <ti{a) with a constant ^1(0;) > 0, and the estimate 

x€ Supp(g) BeB 

r~y 1 (128) 

< U[2(5Y..eB \r}.\ + \\U^m\ for B = {B} 

y 1 otherwise 

which will be enough in order to get the factor in what follows. Using Propo- 
sition 8.3 and Lemma 8.5, we get a uniform upper bound \U'^{B)\ < c\B\ with a 
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constant c > such that c J, for /3 t oo. Hence, in the case B — {B} we have 

\wt{B = {B}) I r ]dAc \ 



2 

xeB BeB BeB 



VdAc 



(129) 



<e-\''h'{2(3+^)\B\E[\U'^{B)\\7joA,] 



Note that the above uniform upper bound on\U'^{B)\is not sufficient to get a sen- 
sible estimate on the conditional expectation. However, a more detailed upper 
bound can be obtained. Without loss of generality, we can assume that B n dAc = 
0, so that the conditioning on tiqac can be omitted. First, assume there is an 
aggregate^ /C^ such that Dom(/C^) = B. Then, Lemma 8.5 gives the estimate 
logZ^ < cyjSl and, since |9/C^| > \B\/2, Proposition 7.5 reads that the probability 
of such an event is bounded by cxp(— y Second, assume there is a family of 
aggregates (/C^Ji (of possibly different orders) such that := UiDom(/C^,) C B. 
Then, any cluster C such that Dom(C) = B has the length \C\ > \B \ D'^\ and 
Proposition 8.3 gives the estimate 

E m<^M-^\B\D''\) 

Dom(C)=S 

Moreover, the probability that — D for a fixed set D is bounded by e"^'^'. 
Note, however, that the above two scenarios are possible only provided that > 
/i, otherwise we only get a contribution from 0-clusters, the sum of which is 

- — IBI 

bounded bye ' ' . All in all, we obtain 

E[\U^iB)\\ < c,|S|e-f l^ll|^|>, +e-4l^l + ^ e'^l^'-^I^V^I 

DCB (130) 

<e-4l^l + l|5|>,^(c7+l)|5|e-^l^l 

provided that f3/lo is large enough. Recall that C5 does not depend on Iq, which 
means that the latter can be adjusted as large as necessary. Using the same argu- 
ment for U~'^{B) and substituting (130) into (129), we get 

E = W)haAc) I e^l^l< 2-^2(2/3 + ^) ^ |i^| e^^+^^l^l 

B3x B3x 
BCdA BCdA K^-'^i-) 



X [e-4l«l + l|5|>,^(c7 + l)|S|e-^l^l] <t'0 



For simplicity, we suppress the subscript n here. 
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which is true for any r' > provided that r and a are chosen sufficiently small 

and /o (and hence l\) sufficiently large. This argument can easily be generalized 
by taking into account all collections B, card(i3) > 1. Hence, the proof of (127) 
is completed by choosing r = under the condition \t\ < to with to = to{(3) 

being small enough. 

12. Concluding remarks and some open questions 

Our result that a typical boundary condition (w.r.t. a symmetric distribution) 
suppresses both mixed and interface states explains why these states are typically 
not observed in experimental situations without a special preparation. To a certain 
extent it justifies the standard interpretation of extremal invariant Gibbs measures 
as pure phases. 

Although this result, which finally solves the question raised in [39], is only 
about the 2-dimensional Ising ferromagnet, and thus seemingly of limited inter- 
est, it is our opinion that the perturbation approach developed in the paper is 
actually very robust (compare [24]). As we have observed at various points in the 
paper, there seems to be no barrier except some technical ones to extend the anal- 
ysis to the Ising model with random boundary conditions in higher dimensions. 
In fact, there might be extensions of our approach into various different direc- 
tions. In particular, both the random distribution of the boundary terms and the 
phase transition itself could lack the plus-minus symmetry, and one might also 
consider a more general Pirogov-Sinai set-up in which the niraiber of extremal 
Gibbs measures could be larger than two. Another possible extension could be 
to finite-range Hopfield-type models, in which periodic or fixed boundary con- 
ditions lack a coherence property with respect to the possible Gibbs measures, 
and thus are expected to behave as random ones [44]. Actually, our result can 
be translated in terms of the Mattis (= single-pattern Hopfield) model with fixed 
boundary conditions, proving the chaotic size-dependence there. 

More generally, in principle the phenomenon of the exclusion of interface states 
for typical boundary conditions might well be of relevance for spin glass models 
of Edwards- Anderson type, which has indeed been one of our main motivations. 
Our result illustrates in a simple way how the Newman-Stein metastate program, 
designed for the models exhibiting the chaotic size-dependence, can be realized. 
The number of states, as well as the number of "physically relevant" states for 
short-range spin glasses has been an issue of contention for a long time. In this 
paper, we have provided a very precise distinction between the set of all Gibbs 
states, the set of all extremal Gibbs measures, and the set of "typically visible" 
ones, without restricting a priori to the states with a particular symmetry. We hope 
the provided criterion might prove useful in a more general context. 
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We mention that the restriction to sparse enough sequences of volumes is essen- 
tial to obtain almost sure results. Actually, for a regular sequence of volumes, we 
expect all mixtures (in dimension three all translation-invariant Gibbs measures) 
to be almost sure limit points, although in a null-recurrent way. This still would 
mean that the metastate would not be affected, and that it would be concentrated 
on the plus and minus measures. See also the discussion in [16]. However, prov- 
ing this conjecture goes beyond the presented technique and remains an open 
question. 
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Appendix A. Cluster models 

In this section we present a variant of the familiar result on the convergence of 
the cluster expansion for polymer models, which proves useful in the cases when 
the summation over polymers becomes difficult because of their high geometrical 
complexity. Such a situation arises, for example, in the applications of the cluster 
expansion to the study of the convergence of high-temperature (Mayer) series in 
lattice models with an infinite-range potential. Since the Mayer expansion tech- 
niques are by no means restricted to the high-temperature regimes (note e.g. its 
application in the RG schemes for low- temperature contour models), the result be- 
low can be applied in a wide class of problems under a perturbation framework. 
In our context, we use the result to provide upper bounds on the weights z/'n of 
n-clusters, see Section 11.2. 

We consider an abstract cluster model defined as follows. Let G = {S, oo) be 
a finite or countable non-oriented graph and call its vertices polymers. Any two 
polymers X ooY are called incompatible, otherwise they are compatible, X ~ F. By 
convention, we add the relations X 'y^ X for all X E S. Any non-empty finite set 
A c S* is called a cluster whenever there exists no decomposition A = Ai U A2 such 
that Ai and A2 are non-empty disjoint sets of polymers and Ai ~ A2, where the 
latter means that X ~ y for all X e Ai and Y e A2. Let V{S) denote the set of all 
finite subsets of S and C{S) denote the set of all clusters. A function g : V{S) ^-^ C 
is called a weight whenever 

i) gm = h 

ii) If Ai ~ A2, then ^(Ai U A2) = ^(Ai) g{A2). 
If the extra condition 

iii) 5'(A) = whenever there is an X G A such that X oo A\{X} 
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holds true, then we obtain the familiar polymer model. In the sequel we do not 
assume Condition iii) to be necessarily true, unless stated otherwise. 

Note a simple duality between the classes of polymer and cluster models: Any 
cluster model over the graph G = {S, ^) is also a polymer model over the graph 
G' — {C{S), ^). The other inclusion is also trivially true. A natural application 
of this duality is to the polymer models with a complicated nature of polymers. 
Such polymers can often be represented as clusters in a new cluster model with 
the polymers being simpler geometric objects. 

To any set A e V{S) we assign the partition function Z{A) by 

ZiA) = J2 9{^) (132) 

AcA 

The map between the functions g and Z is actually a bijection and the last equa- 
tion can be inverted by means of the Mobius inversion formula. In particular, we 
consider the function g'^ : V{S) ^-^ C such that the Mobius conjugated equations 

\ogZ{A) = /(A) /(A) = 5^(-l)l^\^llogZ(A) (133) 

AcA AcA 

hold true for all A e V{S) and A e V{S), respectively. The function g'^ is called a 
cluster weight, the name being justified by the following simple observation: 

Lemma A.l. For any cluster model, g'^{A) — whenever A is not a cluster. 

A familiar result about the polymer model is the exponential decay of the cluster 
weight g^ under the assumption on a sufficient exponential decay of the weight 
g, see [29, 35]. We use the above duality to extend this result to the cluster models, 
formulating a new condition that can often be easily checked in applications. 

Proposition A.2. Let positive functions a,b : S ^-^ R'^ be given such that either of the 
following conditions is satisfied: 

(1) (Polymer model) 
Condition iii) is fulfilled and^ 

(2) {Cluster model) 

There is z : S ^Rf^ satisfying the condition 

sup -i- y e(2"+'')(^);2(F) < 1 (135) 



■^We use the convention § = here. 
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such that \g{A)\ < UxeA z{X)for all A e V{S). 



Then, 

sup , . 
xes a{X) 



E e^-^-'(^^|/(A)| < 1 (136) 



Proof. (1) For the case of the polymer models, see [29] or better [35] for the proof. 
(2) To prove the statement for a cluster model, we represent it as a polymer model 
over the graph {C{S), ^) and make use of the above result. Hence, it is enough to 
show the inequality 

e^yeA(«+^')m|^(A)| <a(X) (137) 



Aec(S) 



for all X e S. Indeed, then one gets 

J2 eS^eA(«+6)(:i^)|^(A)| < J2 a(^) (138) 
Asc(s) yeAo 

for all Ao e C{S) and the statement about the polymer models yields 

A*ec(C(S)) FeAo 

A**-Ag 

where the sum on the LHS is over all clusters incompatible with Aq in the polymer 
model with the set of polymers C{S). Since the weights g'^{A) of the clusters in 
the original cluster model are related to the cluster weights g'^{A*) in the polymer 
model under consideration as 

/(A)= Yl 9^(^*) (140) 

A*: U^,g^,A'=A 

we immediately get 

^ e^^eAK^I^T^^^ll ^ e^AeA*EYeAbiY)\gTi^^*^^^^^x) (141) 

A-toX A*6C(C(S)) 

A*^X 

which is inequality (136). 

Using the notation z{X) := z{X) e«(^)+''(^) and 

for any A e V{S) and X e A, inequality (137) follows from the next two lemmas. 

□ 
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Lemma A.3. The function 3x{A) satisfies the recurrence inequality 

3x{A)<z{X)exp[ J2 MA\{X})] (143) 



YeA\{X} 



Proof For any cluster A we split A \ {X} into connected components, i.e. a family 
of clusters (Aj), and subsequently write: 



3x{A) = z{x) E n n ^(^) 

AeC(A) j YeAi 

ACA\{X} 

oo ^ n 

swE^ E n E 

n=0 yi,...,y„gA\{x} j=i A^cA\{x} ygA,- 

oo ^ 

= E 3r(A\{X})]'^ 

n=0 ' Y*x 

YeA\{X} 

= i(X)exp[ MA\{X])\ 



(144) 



yeA\{x} 



□ 



Lemma A.4. Assume that 



Then 



i(y)e"(^) < a{X) (145) 



J] 3y(^) < a{X) (146) 

Proof. We prove the inequality 

3x(A) < i(X)e«W (147) 

for all A e 'P(<S') and X e A, by induction in the number of polymers in the set A. 
Assuming that this bound is satisfied whenever l^l < n, we can estimate 3x(^) 
for I A| = n by using Lemma A.3, condition (145), and the induction hypothesis as 
follows: 

-bx{A) < i(X)exp[^ 5(r)e'^(^)] < z{X)e<^^ (148) 
y^x 

As the statement is obvious for | A| = 1, the lemma is proven. □ 



KING WITH RANDOM BOUNDARY CONDITION 



51 



Appendix B. Interpolating local limit theorem 

We present here a simple general result that can be useful in the situations where 
a full local limit theorem statement is not available due to the lack of detailed con- 
trol on the dependence among random variables the sum of which is under con- 
sideration. For a detailed explanation of the central and the local limit theorems 
as well as the analysis of characteristic functions in the independent case, see e.g. 
[13]. Here, under only mild assumptions, we prove an asymptotic upper bound 
on the probabilities in a regime that interpolates between the ones of the central 
and the local limit theorem. Namely, we have the following result that is a simple 
generalization of Lemma 5.3 in [16]: 

Proposition B.l. Let {Xn)ne¥s be a sequence of random variables and denote by ipn{t) the 
corresponding characteristic functions, ipnit) = Ee**'^". 7f (yl^jneN/ (5n)neN^?wd (r„)„eN 
are strictly positive sequences of reals satisfying the assumptions 

i) EK^^oo^nH dt\ilJn{t)\ < 27r 



ii) There is k > 1 such that lim„_^oo jjr^ — 

then 



lim 4^ P{a6n <Xn<bSn}<b-a (149) 

n— ►oo On 

for any a < b. 

Remark B.2. Note that: 

(1) Up to a normalization factor, Condition i) of the proposition only requires An to 
be chosen as 

A,,=o(^[r dt\Mt)\y') (150) 

(2) If there is ei such that AnTn < n^^ eventually in n, then Assumption ii) of the 
proposition is satisfied whenever SnTn > rf^ with a constant 82 > 0. 

(3) The choice Sn = An (if available) gives an upper-bound on the probabilities in the 
regime of the central limit theorem. On the other hand, Sn = const corresponds 
to the regime of the local limit theorem. However, for the latter choice it can be 
difficult to check the assumptions, and that is why one has to allow for a sufficient 
scaling of6n, see Part (2) of this remark. 

(4) Much more information about the distribution of the random variables Xn would 
he needed in order to get any lower bounds on the probabilities (except for the 
case Tn = 00 in which a full local limit theorem can be proven). This is a hard 
problem that we do not address here. 

Proof. Let sequences (An), (Tn), (5„) be given such that the assumption of the propo- 
sition is true and take an arbitrary positive function h G C°°(M) for which i) 
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h{x) — for any x ^ and ii) J^^dxh{x) — 1. Using the notation Gn for 

the distribution function of Xn, we consider its 'regularized version' Gn defined 
by the Lebesgue density 



dx 



/oo 
dGn{y)hn{x-y) (151) 
-oo 



where hn{x) := j-h{-^). Obviously, G C°°(]R) and it can be expressed by the 
Fourier integral as follows: 

^^^1- rdte-^'^Mt) rdyeMy) 

ax ZTT J_oo J-oo (152) 



27r/ 



oo 
oo 



^ , dte-''^^|;n{t)h{t5n) 

ZTT 

where h{t) := dx e^^^h{x) and we have used that ^nii) h{tSn} € follow- 
ing from Assumption i) of the proposition and from the bounds |t/^„(t) | , \ h{t) | < 1. 
Moreover, if A; > 1 is such that Assumption ii) holds, then, using the bound 
\h{t)\ < c\t\''^ which is true with some constant c for all t e IR \ {0}, we obtain 
the estimate 

supA„^%^ <^]I^(An r dt\Mt)\ + / dt\h{tSn)\) 

n^oo ^ dx 27rn^ooV J M[-rn,r„] ^ .^ggx 

1 C An 

< 1 + -1 7 TTT^ = 1 

TT k — 1 n^oo d„r„ ^ 

Finally, by using the inequality 

P{aSn <Xn< bSn} = / dGn{y) / dxhn{x - y) 

J aSfi *^ —oo 



/ da; / dGMK{x-y) (154) 

J {a—e)Sn J aSn 

/ dGnix) 

J (a—e)Sn 



we get 



lim --^P{aSn <Xn< b6n} <{b-a + 2£) lim sup A„ — <b-a + 2e (155) 

n^oo On n-+oo ^ dx 

and the proposition follows by taking the limit e ^0. □ 
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